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A quantum anti-ferromagnetic spin-1 model is characterised on a 2D lattice with the following 
requirements; i) The Hamiltonian is made out of nearest neighbour interactions, ii) It is homoge- 
neous, translational and rotational invariant, iii) The ground state is a real singlet state of SU(2) 
(non-chiral). iv) It has a local spin-1 representation. 

Along the way to characterise the system, connections with classical statistical mechanics and 
integrable models are explored. Finally, the relevance of the model in the physics of low dimensional 
anti-ferromagnetic Mott-Hubbard insulators is discussed. 

PACS numbers: 



I. INTRODUCTION 



The description of the low temperature properties of a given physical system is one of the most fundamental problems 
in condensed matter physics. Following the renormalization group ideas the problem of analysing the thermodynamic 
features of a system can sometimes be recast into the characterisation of the infrared states that are compatible with 
the symmetries of the problem. Concepts like the order parameter and the related broken symmetries give us deep 
insights about the properties of the ground state and phase transition between different states. 

So far, universality classes like Fermi-Landau theory appear as paradigms for the classification of quantum systems. 
Nevertheless, as soon as the interactions between the constituents in the system become more and more important, 
the Fermi-Landau picture breaks down and new universality classes or infrared points are needed to study or describe 
the physics of the problem. For instance, quantum states like the BCS [l| or the Laughlin wave functions are 
known to describe the basic features of conventional superconductors or fractional quantum hall systems, respectively. 

In the area of quantum magnetism, the Neel state and the spin wave theory appear as the classic paradigm (see for 
example 0, i, and references therein). In the last decades, a new family of states known as spin liquids [1,0] has 
attracted a lot of attention (see also Q and references therein). One of the reasons for such attention is due to the 
fact that a subclass of these spin liquid states, the valence bond states ( VBS), has been purposed as the right picture 
to describe the low energy properties of Mott-Hubbard anti- ferromagnetic insulators, which have a deep relation with 
the high temperature superconductivity in cuprates and colossal magneto- resistance in manganites . 

From the quantum information community, ideas and tools [l^ . ITsl . [lH have brought new insights and a fresh 
perspective to problems in quantum physics, in general, and highly correlated quantum systems (l6j . in particular. In 
fact, it was pointed out by Preskill[l3| that this interdisciplinary area of research could provide a better understanding 
of questions that appear in many-bod y q uantum entanglement and more refine classifications of different phases that 
emerge in strongly correlated systems [l8|. 

Following this line of research, there has been a strong effort in describing and analysing a class of quantum states, 
known as stabilizer states, in the context of quantum error correction [T9| and quantum computation [20j . Surprisingly, 
this ver y sa me class of states appears in a broad family of spin models that can describe topological quantum theories 
[m, m, H^. As discovered by Kitaev, the understanding of these theories could provide new ideas about hardware 
for self-protected quantum memories and quantum computation [23, . 

In the next sections, we will analyse a class of quantum states with a local tensor description. One of the most 
studied state with this local structure characterises the ground state of the Affleck-Kennedy-Lieb-Tasaki (AKLT) 
model [H, [13], a spin-1 chain with a Heisenberg-like Hamiltonian. The purpose of this work is to extend the AKLT 
construction to higher dimensions and study the physical properties that emerge due to this local tensor structure 
from the quantum state. 

To achieve these goals, several partial results will be presented in this work: 

• An isotropic antiferromagnet spin-1 state is characterised just with symmetry arguments. This step can be 
seen as a generalisation of the AKLT construction to two dimensions where we allow to have multipartite bond 
states. 
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• A parent spin-1 Hamiltonian is defined that has the multipartite valence bond state as a ground state. It is a 
four-body local and real Hamiltonian with explicit SU (2) and parity invariance. 

• The correlations and structure of the multipartite valence bond ground state are analysed using a mapping to 
two dimensional classical statistical models and to one dimensional quantum systems. At this point, we will 
explore the connections of the antiferromagnet state with integrable models. 

• The physical relevance of the model will be motivated comparing it with the interaction and structures that 
appear in usual anti-ferromagnetic Mott-Hubbard insulators. 

The paper is organised as follows: in the next section, we will describe the tensor structure of the ground state of 
the AKLT model. This section will be useful not only to give an overview of the properties of valence bond states in 
one dimension but, at the same time, we will use it to introduce the notation for the rest of the paper. 

The third section contains the main body of this work. In this section, we will define the two dimensional multipartite 
valence bond state and obtain the parent Hamiltonian. We will then use a mapping to a two dimensional statistical 
model to characterise the properties of the ground state. 

Finally, in the last section, we will give a physical motivation of this model. We will summarise some well known as- 
pects and structures, like the hybridisation of the ionic orbitals by covalent mixing, that appears in anti-ferromagnetic 
insulators and that could be described by the valence bond structure of the state we are about to reveal. 



II. ID VALENCE BOND GROUND STATES. 

Many of the features and structures that we will find in this work can be seen as a generalisation of the ground state 
of the AKLT model [26l. [27j. constructed from an antiferromagnet spin Hamiltonian with Heisenberg-like interaction 
between neighbour sites. The basic building block of this ground state is given by a structure, known as valence bond. 
Its inherent quantum nature has brought new perspectives to the field of low dimensional anti-ferromagnetic systems 
extending beyond the model itself. 

Moreover, several features make the quantum phase characterised by the AKLT model, the Haldane phase [1^, [2^ 
important in its own right. First, the Heisenberg model is contained in this phase; this model is the subject of the 
Haldane conjecture which states that systems with half integer and integer spin have completely different behaviours. 
The former are gapless systems while the latter generate a mass gap that cannot be explained with perturbative 
methods and this is a pure quantum effect. Second, this phase is a magnetically disordered phase; it describes 
spin liquids with a hidden topological structure [sot that cannot be represented with the usual two point correlators 
[3i,^,^,iM] (see Fig. [ID. 
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Figure 1; Valence bond ground state. The first line with dark blue spots represents the physical state of the spin chain, while 
the second line represent the implementation of the state with two ancillae systems per site and a maximally entangled state 
between neighbour sites. 

Following the original works [26l. \2^. the AKLT model describes a translational invariant antiferromagnet spin-1 
chain. Its construction is done in two steps: 

1. Every local spin-1 system is constructed from the projection into the triplet subspace of two ancillary spin-| 
subsystems. 

2. In the ancilla picture, two contiguous spin-i of different sites are linked with a singlet state, i.e. a maximally 
entangled state with zero total spin angular momentum (valence bond). 

Denoting the state of a spin-i subsystem by |a) G C^, the second step in the construction of the VBS fix the state 
between neighbouring ancilla spins to 

|o)= E i«)^"/3i/5) = in)-UT), (1) 

{a,P} = {hl} 
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with e|| = — exi = 1 and e|| = = 0. The projection of two spin-i subsystems into the triplet subspace is 

imposed by 

.V2\ + l) a = /3=T 

\i;^^)= {\am + ma))^ \\0) a ^(3 (2) 



V2 



V2|-l) a = /3=i 



In the first part of this equation appears the local state in a tensor notation, in the middle part, the state is described 
in the ancillae spin-i space and in the last part of the equation, we find the same state written in the spin-1 language. 

So that, if we plug in all the local tensor structures to define our state and using a basis in the spin-1 representation 
such that S°'\h) = ze"'"^|c), with {a,6, c} = {a;,y,z} and e°^'^ being the Levi-Civita tensor, the ground state of the 
AKLT model can be described locally as 

m-^)^,= E 7f<7i^)- E (3) 

s={x,y,z} s={x,y,z} 

where are the usual Pauli matrices. This description corresponds to a matrix product state [ssl . [36| . The first 
point to notice is that the characterisation of the state with this particular set of tensors is one of many, due to the 
fact that we can always modify the tensor description with similarity transformations (gauge transformation). So, in 
what follows, we will refer to this particular choice of tensor to describe the triplet SU(2) sector as symmetric gauge 
description (for aesthetic reasons). 

Several features emerge just because of the tensor structure of the state and it can be shown that these characteristics 
remain in the commensurate part of the Haldane phase. 

• Infrared limit.- Although this model is not the infrared fixed point of the phase, it is exponentially close to it and 
contains the long range and topological properties of the Haldane phase. Following Vcrstracte and coauthors, 
a scale transformation can be performed in this class of states. At any point in the scale transformation, the 
local state is written by, 

l*(A)) = ^^a"|0)+^^ E -^1^) (4) 

s={x,y,z} 

where |0) is the singlet state, {|a;), \y), \z)} form the triplet and the correlation length of the state is given by 
^ = i^^i ■ Then, a step in the scale transformation is defined by 

|*(A)) ^ |*(A)), A^^A (5) 
As a consequence of this analysis, we can see how the entropy evolves for any given number of sites L and 



correlation length ^ 



-1 



S{L,A) = 3 log2 log2 . (6) 



Whenever L oo or ^ — > 0, the entropy goes to S* — > logj 4 (see also [13] )■ 

Compositions rules.- We have seen that the tensor structure of the state is due to the fact that we are representing 
the triplet sector of SU{2) as the symmetric projection of two spin-i subsystems. If we relax the projection 
and also allow the singlet to have a non zero amplitude of probability, the VBS is just the decomposition of two 
spin-i into invariant subspaces of SU{2), i.e. ^ (g) ^ = © 1. 

But reading the VBS as a matrix product state, i.e. |a) = J^s p ^ap[s]\s)\P), we find that the composition of a 
spin-1 or spin-0 particle with a spin-i is constrained, i.e. ^ = ^ and ^ ® 1 = i. 

This constraint in the composition of the angular momentum describes an algebraic structure known as quantum 
groups (see for example [38|). In this particular case the quantum group is labelled by SU {2)2- The possible 
relation between the matrix product state construction and quantum groups was already suggested by Sierra 
and Martin-Delgado [1^ 
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• Non local properties.- Finally, we would like to comment about the relation between the hidden topolo gica l 
order of the Haldane phase 13011 and the potential application of this phase in quantum communication [40[. 
Verstraete and coauthors [32l. l33l. [33 | showed that the behaviour of the non-local string order parameter of den 
Nijs and Rommelse [3l| can be understood in terms of the special entanglement properties of the Haldane phase. 
Using the VBS structure of the AKLT model, they showed that it is possible to create a maximally entangled 
pair with local operations and local measurements at any distance. It can be shown that the VBS structure 
does not change in the commensurate part of the Haldane phase and therefore any ground state in the phase 
can be used to obtain a maximally entangled state, at any distance, using local measurements. 



III. TWO DIMENSIONAL ANTI-FERROMAGNETIC SPIN-1 SYSTEM. THE MODEL 



In this section, we will characterise a quantum system in a two dimensional lattice following the valence bond 
construction of AKLT [26l . [27j in one dimension, but in a way complementary to their extension to two dimensions. 
The quantum system is characterised by a Hamiltonian that is composed of nearest neighbour plaquette interactions. 
Geometrically, it is homogeneous, translational and rotational invariant, i.e. invariant under all point group sym- 
metries of the lattice where it is placed. The ground state of the four-body local Hamiltonian is a real singlet state 
of SU{2), i.e. it is a non-chiral state respecting time-reversal and parity symmetries. The local Hilbert space, the 
physical degrees of freedom, are spin-1 particles that live on the links of a square lattice. 




Figure 2: Lattice structure describing the model. At every bond, there is a spin-1 particle. The interactions take place around 
every vertex and involve the four nearest neighbour spins. If the physical degrees of freedom are located at the vertex of a 
square lattice, the system would be represented by a checkerboard lattice. The ground state of the system is characterised by a 
local tensor structure so a convention regarding the indexing of the tensors is needed. The picture shows the convention we are 
going to follow in this work, but it is totally arbitrary. The numbers are used to identify every spin involved in the interaction 
around the vertex. 

The local interaction in this model appears around every vertex in the square lattice and involves up to four spin-1 
particles. The structure of the Hamiltonian, that we are going to present, will be motivated at the end of this work 
where we will connect this model with the structure of the interactions that can appear in usual anti-ferromagnetic 
insulators (see for instance [U, |42| and references therein). In this work, we use a basis of the spin-1 representation 
such that 5^|(t) = ie^'^''|s), with {^,a,s} = {x,y,z} and e^'^'' denoting the Levi-Civita tensor (the completely anti- 
symmetric tensor). The Hamiltonian for this system is described by: 

-f^vertox =gia (^Si ■ S2 + S2 ■ S4 + S3 ■ S4 + Si ■ 63^ + gib (^Si ■ Si + S2 ■ 83^ 

+ 92a ({Si ■ S2f + (5^2 • 5^4)' + (^3 • 5^4)' + (^1 • 83?) + 92b ({Si ■ Sif + {S2 ■ S3) 

+ 93 ({Si ■ Si){S2 ■ S3) - {Si ■ S2){S3 ■ Si) - {Si ■ S3){S2 ■ 54)) 
+ 34a ({Si ■ S2 + S3 ■ Si){Si ■ S3 + S2 ■ Si) + {Si ■ S3 + S2 ■ Si){Si -82 + 83-84 
+ gib I (5*1 • 5*2 + 5*3 ■ 5*4 + >5i • S'3 + >S'2 ■ S'4)(S'i - 84 + 82 - S3 



(7) 



A J2 ({Sl-S2){S2-Si){Si-S3){S3-Si) + {Si-S3){S3-Si){S4-S2){S2-Si) 
perm{1243} 



The first term in the Hamiltonian corresponds to bilinear anti-ferromagnetic interactions between nearest neighbours; 
the second term describes an anti-ferromagnetic interaction between next nearest neighbours at every vertex. For the 
model we arc going to study, its strength is about gib ~ 0.177gia- The third and fourth terms in the Hamiltonian 
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represent similar interactions to the first and second terms but they are biquadratic operators. These kind of terms 
are allowed by symmetry considerations and by the local spin-1 representation of SU(2) operators. Their amplitudes 
are about g2a — 0.312(7ia and 52b — 0.034(7ia. The following terms are several multi-spin exchange interaction needed 
to stabilise the phase we are going to study. Their couplings are of the order of 53 ~ O.OSGgia, .94a — 0.165f/ia, 
gib ^ 0.120.gia and A ~ O.OOSgia. 

The values of the coupling constants are a particular instance of a more general set of Hamiltonians and they just 
give an estimate. We will see in the following sections that the couplings in the Hamiltonian may, in fact, change 
quite drastically without changing the properties of the ground state. This fact gives a hint at the stability of the 
phase we are about to reveal. 



A. The Ground State 

In what follows, we will see that the unique ground state of the Hamiltonian we have just presented corresponds to 
a valence bond state with a four-partite entangled state associated with the vertices of the lattice. The properties of 
the quantum system are inherited by the constraints imposed at the vertex and the dimensions of the local Hilbert 
space placed at every link. As a consequence, the parameters of the tensor structure that characterise the state will 
be completely defined by symmetry arguments. 

We are first looking for a uniform, translational invariant singlet state in the square lattice. Then, we restrict our 
attention to states of the form: 

1. At every link, a local Hilbert space is defined by a projector of two ancillae subsystems in the spin-1/2 repre- 
sentation to the triplet representation of SU{2). 

2. At every vertex, four ancillae subsystems are entangled in a singlet of SU{2). Following [i^ two different states 
are allowed, depending on their chirality (left and right handed). 



S 

o o = 

q k 




+ C 



Figure 3; Representation of the physical Hilbert space defined at every link from tfie projector of two spin-1/2 ancillae 
subsystems and the entangled state at every vertex. In the analysed example, a singlet state is chosen at the vertex and it 
turns out to be a special case of the classical six vertex model. 

Mathematically, both conditions are imposed using a set of tensors. First, at every link, the local quantum states are 
characterised by: 

|^)= (8) 

s={x,y,z} 

with a{s) some complex probability amplitude such that: a(x) = a{y) = a(z) are the probability amplitudes of the 
triplet sector; and a"^ correspond to the usual Pauli matrices. 

At the vertices, four spin-1/2 meet to form a singlet state of SU{2). Decomposing the Hilbert space spanned by 
these spins into invariant subspaces of SU{2), i.e. i®^(8)i(8)i=0+©0_©I©2, we find that the singlet sector 
is not unique but it is two dimensional. This subspace was studied by Wen, Wilczek and Zee. They characterised it 
by an emergent feature of the four spin-1/2 states known as chirality. For example, the singlet sector with positive 
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chirality is explicitly given by'^ : 

|0+) = (l TTii) + I iitt) + Titi) + I itit)] + e-""/=^[| itti) + I 

{ai,b2,C3,d4}G{T,i} 

The singlet state with negative chirality |0_) corresponds to the complex conjugate of |0+). The first point to notice 
is that these two states form a complete complex orthonormal basis for two dimensional singlet sector. Thus, any 
other state with zero total spin angular moment correspond to a superposition of these two states. For example, the 
state made out of two singlets pairs between the first spin and the second one and the third with the fourth spin can 
be rewritten as: 

|0i2034) = 4 E (r:;f,t[o+]-r^fcJ[0-])l«i62C3d4) (lo) 

This equation is a particular instance of more general relations known as skein relations. 

Due to the fact that we are looking for a real SU{2) singlet, the entangled state that defines our system should be 
of the form: 

K\t [0] = ^"^^alit [0+] + e-^*Tllil [0_], (11) 

a superposition of the two chiral sectors with equal weight. 

If we rewrite the state in the vertex in the symmetric gauge, we find that any real singlet state can be described 
by: 

TilillO] = cos0a°^,^<,, + ^ + + (12) 

where </> is the relative angle between the positive and the negative chiral singlet sectors. 

It is straightforward to realise that the multipartite entangled state at every vertex corresponds to a special case of 
the classical six vertex model (see for instance [H, H^l and references therein) , so that the state is decomposed in a 
classical structure (scaffolding) defined at the vertex and the link where we place the triplet degrees of freedom. 

So, in our problem of defining a SU (2) singlet in a square lattice, we found that the vertex of the state is completely 
defined up to a phase. Then, we should impose another physically motivated restriction to fully characterised the 
multipartite entangled state. This constraint comes from the geometric properties of the lattice where the system is 
placed. To make explicit the rotational symmetry that we are going to impose in the state, we will not work with the 
chiral basis; instead we will decompose any possible singlet state of four spin-1/2 by the superposition of two singlets 
pairs between the first spin and the second one and the third with the fourth spin and another two singlets pairs 
between the first spin and the third one and the second with the fourth spin i.e. {IO12) ^ IO34), IO13) (X) IO24)}. This 
basis is a complete basis for the two dimensional singlet sector but it is not orthonormal. Nevertheless, writing the 
state in this basis makes explicit the properties of the state under rotation. 

We arc looking for a rotational invariant state, so the multipartite entangled state at every link should be an equal 
weighted superposition of the singlets IO12) ® IO34) and IO13) (g) IO24), i.e. 

|0) = e*''|0i2) <E> IO34) + e-''^|0i3) <E) IO24) (13) 

Summarising, the state that we are characterising is a real SU{2) singlet that is rotational invariant. Only two 
solutions are possible with the underlying tensor structure we are using. Keeping the symmetric gauge, the two 
solutions are given by: 

1. A set of one dimensional structures along the rows and columns. They correspond to having unentangled AKLT 
chains at every row and at every column. Mathematically the tensor that characterises the state at the vertex 
is given by 



Mm 



(9) 



^ We remind the reader that the physical degrees of freedom are given by spin-1 systems; the ancillary spin-1/2 system can be seen as a 
tool to represent the state. To emphasise the difference we will use curly and normal bra-kets , i.e. \d') and |s) S C'^. 
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2. A multipartite entangled state of the four spin-1/2 ancillae subsystem such that the two states that come from 
the column combine in the triplet sector and then decay along the row. The tensor that represents this state is 



For the first solution, almost everything is known, due to the fact that it corresponds to uncntangled one dimensional 
structures. We will not go on analysing the phase of this state but remind the reader that there are important one 
dimensional structures like stripes states in two dimensional anti-ferromagnetic insulators (see for example [Zsl]). 

In the next sections, we will analyse and study the second possible solution. 



Once, we have characterised an isotropic singlet state in a square lattice, we can look for other models that generalise 
the SU (2) symmetry. For example, in the hexagonal lattice, we can substitute every link with a projector of two 
three level systems that transform as 3 3 into the octet, i.e 3 (E) 3 — > 8. At every vertex, we will find three states 
3(g)3(g)3or3(E)3C3)3. It happens that the SU (3) singlet sector is unique and it is given by the Levi-Civita tensor of 
three indices. So that, the SU(3) singlet state in the hexagonal lattice is completely defined. 

Looking at the geometric constraints in our model, the only requirement that we need in the structure of the lattice 
is that the connectivity of the lattice is of degree four. The square lattice is a particular instance that fulfils this 
requirement, another instance is the Kagome lattice. Nonetheless, in this case, the rotational symmetry on the state 
will come under 7r/3 rotations instead of the 7r/2 rotations of the square lattice. 

Finally, we would like to comment about the associativity properties and consistency condition of the tensor at 
every vertex. Up to now, we were just composing four spin-1/2 systems. If we carry on this composition in steps 
where we join or fuse two subsystems at a time the result can not depend in the way we choose the subsystems to 
fuse. This means that the Hilbert space spanned by the states where we fuse the first and the second with the third 
and fourth spin should be unitary equivalent to the Hilbert space spanned by the states where we fuse the first and 
the third with the second and the fourth spin. The unitary transformation that relates both equivalent spaces only 
depends on the total angular momentum of the states and in our case is given by the classical Racah's symbols (see 
for instance [i^ and references therein): 



where p and q label the total spin angular momentum of the intermediate states in the composition of the four spin-1 /2 
systems (see Fig. H]). 




(15) 



Small remarks and generalisations. 




(16) 



k 





k 



k 




Figure 4; Graphical representation of the unitary equivalent Hilbert space coming from the composition of the different 
subsystems in different order. The unitary that relates the Hilbert spaces is given by the classical Racah's symbols. 
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B. Local spin Hamiltonian 

The quantum system that we have just defined in the last section is characterised by the state: = 
X]s={a: y z} f I*) ^t cvcry link and the tensor: r[0] = (cr^ ® + ® + iSi a^) at every vertex. Around 
every vertex, this tensor structure spanned a sixteen dimensional subspace JC from the eighty one dimensional Hilbert 
space of four spin-1. Any positive definite operator that involve the four spin-1 states around every vertex, and whose 
kernel is exactly the subspace /C, is a possible Hamiltonian, for which the valence bond state is the ground state. 
As wc have seen, symmetry and geometric constraints fully specify the valence bond state. The freedom in choosing 
the Hamiltonian gives an idea of how stable this kind of valence bond phases are under any perturbation. In fact, 
the Hamiltonian need not respect the SU{2) symmetry and still the ground state will be given by the valence bond 
construction. 

Due to the interactions that usually appear in anti-ferromagnetic insulators, we find that the Hamiltonian given by 



i^vcrtcx —gia (^Si ■ S2 + S2 ■ Si + S3 ■ S4 + Si ■ S3) + gib (^Si ■ Si + S2 ■ Si 

+ 92a [{Si ■ 82)^ + (S2 ■ Sif + {S3 ■ Sif + {Si ■ S3f) + 526 ((5^1 • 5^4)' + (^2 ' S3) 

+ 93 ({Si ■ Si){S2 ■ S3) - {Si ■ S2){S3 ■ Si) - {Si ■ S3){S2 ■ Si)) 
+ 54a ^("^i ■ S2 + S3 ■ Si){Si ■ S3 + S2 ■ Si) + {Si ■ S3 + S2 ■ Si){Si ■ S2 + S3 ■ Si) 
+ 9ib I {Si • 5*2 + 5*3 ■ 5*4 + S*! • S3 + S2 ■ S'4)(S'i ■ 5*4 + 5*2 • 5*: 



(17) 



A ((5l-52)(52-54)(54-53)(S'3-5^l) + (5i-53)(53-54)(54-52)(52-^i) 

porm{1243} 



can be of some relevance in the studies of this kind of materials |4lL |42| . 

Using the parameter A > as a perturbation, the operator _ffvcrtox(A) is positive definite with a kernel (ground 
state) given by the valence bond state, if the other parameters are given by 

1 + 86A 2, 446A-1 3 + 167A 3 + 734A, . 
5ia = l; .916 = ; .g2a = y (1 + 44A) ; 526 = ^ ; 53 = 20A; gia = — ; 546 = (18) 

Moreover, the gap to the first exited state of the local Hamiltonian in one vertex is given by: gap = 360 A. So that, 
when the local gap is of order one, the parameters gi take the values displayed in the introduction. 

Following the same steps as in the AKLT construction [2^, [l^, [13], we can show that the valence bond state is, 
in fact, a unique ground state of -ffvoitex- Nevertheless, there is a sufficient condition to prove the uniqueness of the 
ground state [48[ . The condition states that if the mapping that goes from the Hilbert space spanned by the ancillae 
subsystem in the boundary to the Hilbert space of the physical spin-1 degrees of freedom in the bulk is injective then 
the ground state is unique. 

Once we have specified the state by the tensors A[s\ and r[0] is straightforward to check the injectivity of the map 
and then, the uniqueness of the ground state. 



C. Ground state properties and correlations. Expectation values 



We now want to characterise the state we have just defined and to get its correlation values. The value of any 
expectation value is obtained via a mapping of the 2D quantum state to a 2D classical statistical model and from 
there to a Id quantum mechanical problem using a transfer matrix defined from the 2D quantum state. The long 
distance properties of the state happens to be described by an integrable statistical model. 



1. 2D Statistical Model. - 



In order to obtain any expectation value given a valence bond state, a vertex matrix R is defined from the tensors 
at the links and at the vertices. So that, from these kind of valence bond models, two different classical structures 
appear in the vacuum-to-vacuum expectation value: 
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• At every bond, a transfer matrix is placed, defined by: E — J2s={x y z} kC^)!^ {'^^T ® '^^ 

• At every vertex, a vertex matrix is placed, defined by: 1/ = F* ® F 
and the R matrix is given by: 



(19) 



{jl J2 ,j3 J4} 



For instance, the vacuum to vacuum expectation value in the valence bond state correspond to the calculation of 
the classical partition function of the vertex model: 



(20) 



all configurations lattice 

where one sums over all possible configurations compatible with the weights i?. 
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Figure 5: Vertex model in a classical model. The top part represents the vertex 7?^^, the Boltzmann weight of the vertex model. 
The bottom part represents the row-to-row transfer matrix for the vertex model. 



So, given any valence bond state model described by local projectors in a two dimensional lattice, it is always 
possible to build a classical vertex model defined by a matrix R that can be interpreted as Boltzmann weights of a 
classical model. Then, we have a mapping from a class of quantum states to classical models in a lattice with the 
same dimensions (see also [i^ [s^ [5l|, Is^ ) . 

Once, we have realised about this equivalence between 2D valence bond model and classical statistical models, we 
can use well-known tools from statistical mechanics to obtain information about the quantum state. For example, 
the sufficient conditions for the integrability of a classical model i.e. the star-triangle or Yang-Baxter equations 
[38, , 44, . .531] , will become sufficient conditions to obtain information about the thermodynamic limit from the quantum 
valence bond state. Then, a possible strategy is to study the possible valence bond structures that give naturally 
integrable models. 

Instead of going on with the two dimensional classical analogy, we will use another well-known mapping, the 
equivalence between D-dimensional classical statistical models and d = D — 1 dimensional quantum mechanical 
problems. This mapping is at the core of path integral 5J| and quantum Montecarlo methods [55|. 



2. Id Quantum System. - 

In quantum mechanical problems, time evolution is carried out by the row-to-row transfer matrix T, with matrix 
elements Tab given by 

So, if we understand this time evolution operator as the exponential of a quantum Hamiltonian Hid = — logT then the 
classical partition function Z2D is just given by the trace of the density matrix of a Gibbs ensemble Z2D = Tr I e^^-^" ) . 
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If we are able to diagonalise the one dimensional quantum Hamiltonian Hid, we will be able to extract information 
about the 2D valence bond state. For example, consider the normalisation of the valence bond state, 

im = 22D = Tr(e-^^") = Tr |/i)(A^)^(/x| j (22) 

where | /i e N} is a complete set of eigenvectors of the one dimensional quantum Hamiltonian and are the 
eigenvalues of the transfer matrix T. Then, in the thermodynamic limit iV >> 1, the value of the highest eigenvalue 
is fixed by normalisation to Ai = 1. 

The other eigenvalues give information about the correlation in the system. For instance, the two point function in 
a 2D lattice is given by 

{0,0s) = im (EI'")(^m)"''"'(m|) Og|l) 

(23) 

= E(1|O.Hm)(m|OoH1) (^^j e~, 
where — , defines the correlation length and Og is the classical matrix similar to the classical transfer matrix 



ogJA^I 

T but, in the x-link, we substitute the matrix 



by the matrix 



E= K^)l' i'^r®'^' (24) 

s={x,y,z} 



a{sya{^l) {a')*®a'' {s\0\^i). (25) 

{s,fj.} = {x,y,z} 



For the valence bond model that we are analysing, the local dimension of the classical indexes in the transfer matrix 
is four dimensional. This vector space can be understood as the direct sum of the singlet sector |0) and the triplet 
sector {|a;), \y), \z)}. To built the one dimensional quantum Hamiltonian from the normalisation of the 2D quantum 
state, we will use the spin operators in a four dimensional representation, such that, 

[J\JP]=^e''P^J1^ jfc|0) = o, Ap)^te'''"'\q), {k,p,q} G {x,y, z}. (26) 

and the creation and annihilation operators, 

bl\q)^dgo\k), 5fc|p) = 5fcp|0) {k,p}e{x,y,z} (27) 

With these definitions, the one dimensional quantum Hamiltonian can be recast into, 



+"5 Y {^Lkbn+i.k + fe«,fc6i+i,fc) + ae Y {^i,kbi+i,k + bn,kbn+i.k) + ar 

k—{x,y,z} k—{x,y,z} 



(28) 



with 



logs 31og( sT^^j log (469176871410737- 130126251314344 Vl3) - 20 log 3 



r 1/497+136/13 \ , rT ( 497+136yi3 \ , 

&iog(^ 81 j l0g3^ ^510g(^ j lQg3^ ^51og(^ 



497- 136V 13 



"2 = /— 03 = „ 1 aA = 



log 3 



12\/13 6 8V13 4 12V13 



(29) 



The first two terms in this one dimensional quantum Hamiltonian correspond to a bilinear-biquadratic anti- 
ferromagnetic interactions; the third, fourth and last term are diagonal terms that lift the onsite energy of the 
different states; the fifth and sixth are the hopping and pairing terms, respectively. 
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As we have seen, the spectrum of this one dimensional quantum Hamiltonian gives the correlation function of the 
2D valence bond state. To solve it, we are going to use a numerical renormalization group method. For a related 
work where the renormalisation transformation is applied on the 2D lattice see [s^ . 

But before going into the numerics, we would like to analyse from the structure of the 2D quantum state the 
extreme cases that we could obtain in our numerical study. 

The 2D quantum state, as we have defined, has two structures: at every vertex, there is a tensor r[0] = 
(cr^ cr^ + (T^ (81 cr^ + (T^ (g) (T^) that is equivalent to a particular instance of the classical six vertex model; at 

every link, there is quantum state X){s} ^I^lk) ~ Ss={£c y z} '^^^^\■^)^ that corresponds to a combination of states 
in the triplet sector. To analyses the extreme cases, we introduce a parameter A dependence in the state at every 
link, such that 



|V'(A))=a(0)aO|0)+ a(5)a-^ \s) = ^^i±^a°|0) + z^^i^ \')^ (30) 

s={x,y,z} s={x,y,z} 

where the state |0) represents the singlet state and \z)} the triplet sector. The scalars a(0), {a(s)} give the 

amplitude of probability in the singlet and triplet sectors, respectively. 

When A — — i, we recover the spin-1 model we are dealing with. If A = 1 the whole 2D valence bond state is 
decomposed in an uncorrelated classical tensor at every vertex and a product of singlets states at every link. And if 
A = 0, the valence bond state is decomposed in a correlated classical tensor at every vertex and an equal superposition 
of the singlet and triplet sectors at every link. In this last case, although the classical and quantum structures are 
correlated and there is non-trivial entanglement, the whole quantum state has exactly zero correlation length. This 
feature is due to the fact that the transfer matrix that appears at every link when we calculate the vacuum to vacuum 
expectation value is given by, 

E= E K)*®'^^- (31) 

M={0,a:,y,z} 

If we diagonalize this matrix we find one eigenvalue cq = 1 and three eigenvalues e^- = = 62 = A. This means, 
when A = 0, the transfer matrix is a one dimensional projector and any correlation length is zero (which does not 
mean that the state is trivial, since there are still non trivial quantum correlations, i.e. entanglement). Then, in 
principle, we could think to use the parameter A to apply perturbation theory. The numerical results, that we will 
discuss in the next lines, show that when the A insertion in the links is not a trivial one, the long range behaviour 
of the quantum system corresponds to a conformal invariant theory, i.e. an instance of an integrable model with an 
infinite number of conserved quantities and infinite correlation length. This result should be linked with a theorem 
proved by Hastings [s^l where he shows that if a state has an algebraic decay of the correlation function then the 
Hamiltonian for which this state is the ground state is gapless. 

In the numerical analysis, we use a renormalization group method (CORE) [58l.[59t that combines the decimation of 
the short distance or high energy degrees of freedom with a matrix product state ansatz (ssl . [36j for the ground state. 
The choice of the ansatz is due to the fact that the structure of the one dimensional Hamiltonian Hid is equivalent 
to the interactions that appear in spin chains in the Haldane phase [28l. [2^. There, we know, that the infrared limit 
of the phase is given by the AKLT model, that is a special case of a matrix product state; and the phase transition 
is described in the continuum limit by a Wcss-Zumino-Novikov-Wittcn (WZNW) model on the SU{2) group at the 
level fc = 2 [63,[6ll,[6l. 

Due to theoretical analysis and scaling arguments, we know that, in this conformal model, the minimum gap behaves 
as A(A) ~ N^2^^^\ with 9 = —1, as we modify the infrared cut-off N = 2'"^^'''. FiglHlshows the plots for the behaviour 
of the minimum gap. The slopes in the fitted curves show the scaling of A (A) as the infrared cut-off is modified with 
9 ~ -0.99(4). The plots correspond to values of A = {1/100, 1/30, 1/10, 1/3}. 

In the second plot, Figl?! we check how change the starting point in the minimum gap as we change the value of 
A € [1/100, 1/3]. The fit curve is just to guide the eye but follows a functional form /(A) ~ ^ -|- a2A^= | logs ^1^' 
for some constant values {ai,f3j} given in the footnote of FiglT] 

Before closing this section, we would like to point out the connections between some effective (gauge) theories and 
the features of the model we have just revealed. The structure of this anti-fcrromagnet is given by a vertex tensor 
r[0] and a link state IV^A)) = '^^^^^ a°\0) + i^^^ J2s={x.y.z} k>- When the insertion of the parameter A = 0, 
all the correlations in the system are zero, i.e. the Hilbert space of the theory at the boundary is one dimensional, 
the local Hamiltonian and dynamics can be expressed by a set of commuting projectors, i.e. local constraints, and 
nevertheless it can be said that the theory is by no means trivial. As soon as the insertion of the parameter A is 
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LOG [2, A] 




Figure 6: Plots of the logarithm in base two of the gap from the one dimensional quantum Hamiltonian as short distance degrees 
of freedom are integrated. By scale, we mean the number of steps in the renormalization, where we take groups of two contiguous 
sites every step. The points correspond to the numerical data, from top to bottom, they belong to A = {1/100, 1/30, 1/10, 1/3}, 
respectively. The solid line corresponds to the fitted curve and the dashed line to a 95% confidence interval for the predicted 
responses given a linear fit. The slope of the curve corresponds to a critical exponent 6 ~ —0.99(4). 



05 0.1 0.15 0.2 0.25 0.3 



Figure 7: Plot of the function /(A). The points correspond to the numerical data in the interval A € [1/100, 1/3], the line is 

0.' 



just to guide the eye following the function /(A) ~ + 1.1A" '^| log2 A|^'^. For more details, look in the main text. 



different from zero the theory that appears at the boundary of the two dimensional system is given by a very specific 
intcgrablc model. 

In 1989, Witten [gI, [g^, [1^ noted that the boundary behaviour of a class of quantum field theories known as 
Chern-Simons theories in 2 + 1 dimensions are closely tied to a class of conformal field theories called Wess-Zumino- 
Novikov- Witten (WZNW) models in 1 + 1 [gI, [G^I- Among many other things, Witten showed that the Hilbert space 
at the boundary in these 2 + 1 theories is one dimensional if there is no inclusion of Wilson lines. As soon as this 
inclusion is not zero, it starts to appears the whole Hilbert space structure of the conformal theories (see also [gl, iGOj 
and reference therein). 

If the A insertions in our microscopic model are understood as inclusions of Wilson lines, the system that we have 
just defined and characterised can be understood as a microscopic model of a Chern-Simons theory in the bulk with 
a SU{2)2 WZNW model living on the edge. 

It might be interesting to go deeper into the analogies that appear between the microscopic model, we have studied, 
and Chern-Simons theories. For instance, features like the degeneracy of the Hilbert space depending on the topology 
where is defined should, in principle, be understood from the microscopic anti-ferromagnet. 

IV. ANTI-FERROMAGNETIC MOTT-HUBBARD INSULATORS 

The purpose of this section is not to make any claims about the understanding of the physics of anti-ferromagnetic 
Mott-Hubbard insulators. It merely shows that the ground state structure of these systems could be described by a 
valence bond states in two dimension. A complete comparison with experiments is needed (see for example (70j and 
references therein). 

Most anti-ferromagnetic Mott-Hubbard insulators arc transition metal compounds, in which the d-orbital cations 
are separated by large anions [4l|, An example of these anti-ferromagnets is given by ceramic cuprate systems, 
basic building blocks of high-temperature superconductors. These cuprates are build of layers of Cu02 where the 
copper cations are placed at the vertices of a square lattice separated by oxygen ions. Electrons in the last shell of the 
copper has a partially filled Sd^ configuration, while the oxygen sites are completely filled. Direct hopping between 
(i-orbitals seems rather unlikely. 



Anderson 71| showed that the concept of kinetic exchange is enhanced by mixing between cation and anion orbitals, 
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explaining the large ordering temperatures of these materials. The anion- mediated exchange described in what follows 
is known as super-exchange. To understand how it works, we will consider a simple model that contains the essential 
features to explain the physics in this system. 

The model is made out of two cations with one rf-electron each and the intervening anion is an anion. In the 
ionic picture, i.e. when there is no overlap between the orbital wave- functions of the ions, the oxygen has a filled 
p-shell. As soon as the d and p orbitals overlaps, covalent mixing allows the p-electrons to partially reoccupy the 
cations. Following [t^ I the problem is treated from the point of view of the p-electrons, for the two cases when the 
(i-spins are either in the triplet (parallel) or in the singlet (anti-parallel) configuration. 

The on-site energies of the p and d orbitals are given by Ep ~ {ay\H\ap) and Ed = (adlHlad), with a = {t, i}- The 
atomic levels are understood to arise from the solution of the many-electron problem at the ionic level. We assume 
that Ep < Ed- The covalent mixing amplitude of the p-d orbitals is given from perturbation theory by A ~ '^^e~§~- 



t 

cc ^ 

t 



Figure 8; Lattice structure of usual antiferromagnet insulator. The anion is place at the link of the lattice with a p-orbital 
configuration. The cation is place at the vertex of the lattice in a rf-orbital configuration. Super-exchange between two cations 
takes place via an intermediate anion. In the ionic picture, the anion has the outermost shell filled. Covalent mixing allows the 
anion p electrons to partially reoccupy the cation d orbitals. The mixing depends on the spin configuration. 

If both d-spins are parallel, e.g. | tt); the |-spin p-electron can extend to the cations on its left and right side, 
while the t-spin p-electron must stay at the anion due to the Pauli exclusion principle. 



Ip) covalent-mixing 



\ip) + MU) + MU) 



\p) Pauli principle | tp) 



VI + 2A2 (32) 



The energy of the two p-electrons is given by 



^TT ^ TrWt^^^f I + I + ^^^'^^ I) 1^1 (I W + M id,) + A| idn))] 

^ 1 + 2A^ ~ 

If the (i-electron spin configurations is | jt), the i-spin p-electron can extend to the right cation and the |-spin 
p-electron to the left one, 



I Ip) covalent-mixing 

I tp) covalent-mixing 
The total p-electron energy contribution is estimated as 



I ip) + A| irfj 

x/TTA^ (34) 



+ (35) 
2A 
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The supcr-cxchangc splitting is approximately 

^""^^T-(1 + 2A2)(1 + A2)- 

The parallel configuration costs more energy, thus it is more probable to find the system in the singlet configuration 
than the triplet one at every link. As soon as the orbital wave functions of the different ions start to overlap, a 
splitting between the singlet and triplet sector appears. 

In the valence bond state, we have found that the structure of the state at the vertex of the 2D lattice is completely 
defined by global symmetries. The vertex tensor 

r[</.] =cos0aV + ^ J2 (37) 

defines locally a real singlet state of SU{2). The structure of the state at the link depends on a parameter A that 
describes the splitting in the amplitude of probability of finding the system in a triplet or singlet configuration; when 
A = both sectors have the same amplitude of probability and any correlation in the system is zero. As soon as 
A 7^ there is a singlet-triplet splitting and the 2D anti-ferromagnet starts to develop correlations. 
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